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Eva-Maria Graefe, Maria Graney, and Alexander Rush 
Department of Mathematics, Imperial College London, London, SW7 2AZ, United Kingdom 

We consider a simple quantum model of atom-molecule conversion where bosonic atoms can 
combine into diatomic molecules and vice versa. The many-particle system can be expressed in 
terms of the generators a deformed SU{2) algebra, and the mean-field dynamics takes place on 
a deformed version of the Bloch sphere, a teardrop shaped surface with a cusp singularity. We 
analyse the mean-field and many-particle correspondence, which shows typical features of quantum- 
classical correspondence. We demonstrate that semiclassical methods can be employed to recover 
full many-particle features from the mean-field description in cold atom systems with atom-molecule 
conversion, and derive an analytic expression for the many-particle density of states in the limit of 
large particle numbers. 
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I. INTRODUCTION 

The experimental progress in confining and manipulat¬ 
ing cold atoms and Bose-Einstein Condensates (BECs) 
offers a unique opportunity to investigate the quantum 
properties of interacting many-particle systems. For 
most realistic setups, however, a theoretical full many- 
particle description is beyond the current state of the 
art. Most commonly the mean-field approximation is 
applied, resulting in a description of the many-particle 
system by an effective single particle wave function. The 
time evolution is in this description governed by a nonlin¬ 
ear Schrodinger equation. This approximation is closely 
related to the classical limit of single particle quantum 
systems, where the particle number plays the role of li~^. 
Recently there have been several studies demonstrating 
that this analogy can be used to apply semiclassical tech¬ 
niques to recover full many-particle features from the 
mean-field description alone m- 

Over the last decade there has been considerable in¬ 
terest in atom-molecule conversion in cold atoms and 
BECs [TUIT]. Theoretically these systems are closely re¬ 
lated to deformed SU{M) algebras [THj. Neglecting quan¬ 
tum fluctuations in the many-particle dynamics leads to 
mean-field approximations defined on phase spaces with 
non-standard geometries, and new interesting features in 
the many-particle and mean-field correspondence. Here 
we address the question whether semiclassical methods 
might be adapted to describe the full many-particle 
behaviour of atom-molecule conversion systems on the 
grounds of their mean-field approximation. 

We focus on the simplest model of an atom-molecule 
conversion system, similar to the one considered in 
[7], consisting of non-interacting atoms and diatomic 
molecules each of which can populate only one mode. In¬ 
troducing the atomic and molecular creation and annihi¬ 
lation operators, a’^, a and b\ b, respectively, this system 
can be described by a Hamiltonian of the form 

H = CaCi^a + etPb-\ - -^{a^a^b + aaP), ( 1 ) 

2v N 

where ea,& is the energy of the atomic or molecular mode. 


and V describes the conversion strength between atoms 
and molecules. The total number of atoms N = a^a+2Kb 
is a constant of motion. For a fixed value of the particle 
number N the system lives on an [y] ^ + 1 dimensional 
Hilbert space. For simplicity, we confine the discussions 
to even particle numbers in what follows. 

We begin with a review of the many-particle and 
mean-field descriptions, and their correspondence. Then 
we introduce a semiclassical quantisation condition and 
demonstrate that the many-particle spectrum can be ac¬ 
curately recovered from the mean-field dynamics, and 
derive an analytic expression for the many-particle den¬ 
sity of states in the semiclassical limit of large particle 
numbers. We also present results for the semiclassical 
many-particle eigenstates. We end with a summary and 
outlook. 


II. THE MANY-PARTICLE SYSTEM 


Similar to the two-mode Bose-Hubbard model m , we 
can apply a Schwinger-type transformation to introduce 
the operators 
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( 2 ) 


The Hamiltonian can then be expressed as 


H = eK, + vk,, (3) 

where we have shifted the zero energy and introduced the 
parameter e = 2ea — e-b- 

The physical meaning of the operators Kj is similar to 
the two-mode Bose-Hubbard case. That is, Kz measures 
the population imbalance between molecular and atomic 
mode, and Kx^y measure their phase relation. It is further 
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FIG. 1. (Color online) Eigenvalues of (and Ky) for N = 
50. 


convenient to introduce the operators 

K±=K,±iky, (4) 

that create a molecule out of two atoms and vice versa. 
The operators ([^ and ([^ are related to a nonlinear de¬ 
formation of an SU{2) algebra [UllinilH]- That is, they 
fulfil commutation relations of the form 

[k,,k±] = ±k±, (5) 

just as for SU{2), and 

[k+,k_] = F{k,,N), (6) 

where F(kz,N) is a polynomial in k^ and N. Specifi¬ 
cally we have 

Fikz,N) = -^-^(^N + Akz)[N-12kz). (7) 

In terms of Ky, and the commutation relations 
read 


[kz,ky]=iky, ( 8 ) 

[ky,kz\=iky„ (9) 

[k.,ky\ = ]^F{kz,N). (10) 


The total particle number N commutes with Kx, Ky, and 
kz- In the case of a Bose-Hubbard dimer the conserved 
particle number reflects the conservation of the total an¬ 
gular momentum. For deformed SU{2) algebras this is 
replaced by a less trivial conservation law. We can find 
the conserved quantity using the approach in [T5] as 


C = K_K+ + —K 
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FIG. 2. (Color online) Many-particle spectrum in dependence 
on e for v = 1 and = 10 (top) and TV = 50 (bottom) 
particles. 


Evaluating (C) in any eigenstate of kz yields the conser¬ 
vation law 

2(7V,) (Nkz) (N^kz) 

(Kx) + (Ky) - - — + 

(Nkf) Ajk^z) , N\ N 
N N 16 4 ■ ^ 

We shall see later that this corresponds to a deformation 
of the familiar Bloch sphere of two-level systems to a 
teardrop shape on which the mean-held dynamics take 
place. 

From the commutation relation ([^ it follows that 7V+ 
and K- are the usual lowering and raising operators for 
Kz, from which one can deduce that the spectrum of Kz 
is equidistant [18] . For a given even particle number TV 
the eigenvalues of kz run in integer steps from — to 

N 

4 ■ 

While due to symmetry the operators kx and ky are 
isospectral, their spectrum differs from that of Kz, and in 
particular, is not equidistant. The spectrum is depicted 
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for the example TV = 50 in figure It can be seen that 
the eigenvalues are symmetric around zero, and approx¬ 
imately equidistant at the boundaries of the spectrum, 
while they are closer together around the centre. The 
eigenvalues can be obtained analytically using a Bethe 
ansatz approach [12]. 

Figure depicts the eigenvalues of the many-particle 
Hamiltonian (§ as a function of e for u = 1, and two 
different particle numbers. In comparison to the familiar 
Landau-Zener behaviour of a many-particle two-level sys¬ 
tem without atom-molecule conversion, the simple sym¬ 
metry with respect to e is lost, and we observe a cluster 
of narrow avoided crossings for intermediate values of e, 
stretching from the upper end of the spectrum for neg¬ 
ative values of e to the lower for positive values. These 
are the reflection of the accumulation of eigenvalues in 
the spectrum of Kx in figure For large values of e 
the many-particle spectrum is dominated by the equidis¬ 
tant spectrum of Kz- The narrow avoided crossings lead 
to quasi-stationary states related to unstable stationary 
states in the mean-field description. These are related to 
the molecular mode, where all particles are paired up in 
diatomic molecules. 

The many-particle dynamics can be straight-forwardly 
obtained by integration of the Schrodinger equation on 
the (y + l)-dimensional Hilbert space. It is nevertheless 
instructive to study the Heisenberg equations of motion 
for the dynamical variables Kx^y^z, given by 


f^Kx= -eKy 

^^ky= ekx + '^ + ^ ( tv 2 - skzN - 48k^) (13) 
iK.= vKy 


Due to the nonlinearity of the commutator of Kx and 
Ky this is not a closed set of equations if u ^ 0, that is, 
the right hand side contains dynamical variables such as 
(iV^) whose dynamics is not determined by the dynam¬ 
ical equations (13). In the trivial case v = 0, where the 
atomic and molecular mode decouple, we have a closed 
set of equations, describing rigid rotations around the z- 
axis. In the general case, taking the expectation value 
and neglecting quantum fluctuations, i.e., approximating 
expectation values of products with products of expecta¬ 
tion values, yields the mean-field approximation we shall 
discuss in the following section. It is interesting to note 
that the same dynamical equations in terms of generators 
of a deformed SU{2) algebra also appear in the context 
of fermionic atom-molecule conversion uniiia. 


replacing expectation values of products with products of 
expectation values in the dynamics of the operators Kj 
in equation (13). Introducing the mean-field variables 

Sj = ri{kj), with 77=(^-|-l) ^—>-0, the mean-field 
dynamical equations become 


Sx - 

Sy= eSx -|- — (l — 4 :Sz — ISSj) (14) 

Sz= VSy. 

These can be formulated as canonical Hamiltonian dy¬ 
namics with the classical Hamiltonian function 


H = ri{H) = esz+vsx, (15) 

and the Poisson brackets 

{Sa;) Sy} = — [l — 4Sz — 

= -Sa; (16) 

which directly follow from the many particle commuta¬ 
tors with the identification 


{sj,.Sfc}= lim ir]{[Kj,Kk]). 

N^OO 


(17) 


The dynamics (14) is confined to a two-dimensional 
surface given by the constraint 


sl+4 = \(^~ + 2sz)^ =: k{sz), 


(18) 


with 


e - 


1 ii 


This constraint also follows from the 


many-particle conserved quantity (12) in the mean-field 
limit. The resulting surface is depicted in figure It 
has a characteristic inverted teardrop shape, with a tip 
at Sz = — Note that similar shapes are known as 
Kummer shapes in the classical description of coupled 
oscillators of different frequencies |2^25j . 

From the dynamical equations it follows that all 
fixed points are located at Sy = 0, and have to fulfil 


eSx — ^{^+Sz){2 12sz)- 


(19) 


Using the constraint (18) this yields a polynomial in the 
Sz component of the fixed points 


(o + ■s^)^ 


— k — — 2e4sz + 


Qv'^sl =0. (20) 


III. MEAN-FIELD APPROXIMATION 

In this section we review the derivation of the mean- 
field dynamics, and summarise its most important fea¬ 
tures [HHni' Here we formulate the approximation in a 
way that is most natural from the perspective of a semi- 
classical limit. The mean-field approximation can be ob¬ 
tained in the limit of large Hilbert space dimension by 


It can be verified that for each of the solutions Sz € 
[—of (20) there is one corresponding value fo Sx 
such that the dynamics is stationary. That is, there is al¬ 
ways a fixed point at the tip of the teardrop, at = — ^. 
Depending on the values of v and e there can be one or 
two further fixed points. For values of e that are smaller 
than a critical value |ecrit| = there are three fixed 

points. A stability analysis reveals that the fixed point 
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FIG. 3. (Color online) Mean-field dynamics on the deformed 
Bloch sphere for v = 1 and e = 0, 1, 2 (from top to bottom). 
The right panel shows the dynamics of Sz for selected initial 
conditions. 


at the tip of the teardrop is a saddle point in this case, 
while the other two fixed points are elliptic. At the criti¬ 
cal value Iccritl = V^\v\ one of the other two fixed points 
collides with the fixed point at the tip of the teardrop and 
then moves to unphysical values Sz < — |. The two fixed 
points interchange their stability in this transcritical bi¬ 
furcation. That is, for values of e larger than the critical 
value there are two elliptic fixed points, one of which is 
located at the tip, the other elsewehere on the teardrop. 
Examples of the dynamics for v = 1 and different values 
of e are shown in figure]^ In the top two figures the fixed 
point at the tip of the teardrop is a saddle point, which 
is approached asymptotically by trajectories starting in 
a particular region on the teardrop, that shrinks as e gets 
closer to the critical value. In the lowest figure the value 
of e is super critical, and we observe the typical oscilla¬ 
tions for trajectories in the neighbourhood of the elliptic 
fixed point at the tip of the teardrop. 

We can introduce a set of canonical variables p and q 
via the transformation 

Sx = r{p) cos(g) ( 21 ) 

Sy = r{p) sin(g) ( 22 ) 

Sz=p, (23) 

with r{p) = 2pKi+2p) ^ p g [“I; 5 ] and q G 
[0, 27r]. It is straight forward to verify that this indeed 


recovers the Poisson brackets (16) with the standard def¬ 
inition 


dAdB dAdB 


(24) 


In terms of p and q the dynamics is then given by the 
canonical equations q = ^ and p = with the 


Hamiltonian function (15) expressed in terms of p and 

q- 


H = ep + - 2p)(l -h 2p)2 cos(g). (25) 


The mean-field dynamics can of course be expressed in 
terms of a nonlinear Schrodinger equation for an effec¬ 
tive single-particle wave function. In the case of atom- 
molecule conversion there are two natural ways to define 
a mean-field wave function. The first, most commonly 
used, arises via the usual identification of the components 
of the mean-field wave function with the probability am¬ 
plitudes to be in one of the states, in this case to be 
in the atomic or molecular state. That is, we make the 
identification —?► and —>■ The 

resulting mean-field wave function is then normalised as 
IV'aP + 2|V'bP = 2, where |^ap/2 is the probability to 
find the mean-field system in the atomic mode and |'i/' 6 p 
the probability to find it in the molecular mode. The 
mean-field dynamics is then governed by the Hamilto¬ 
nian dynamics 


ihtpj 


dH 

w 


(26) 


with the Hamiltonian function H expressed in terms of 
the V'j as 


H = \ {Ha? - SiV'bP) + ^ , (27) 


that is 



With the identification 

Sx = ^ {i’f'tpb + V’bV'a) : 

Sy = (V'f V'b - V’bV'a) : (29) 

Sz = ^ {\^a? - “^Hb?) ■ 

this yields the dynamical equations ( |l4| ) as expected. 
Alternatively, we can replace —>■ ^ , and 

^ xi*'’- Then we have \xa\ + 2|xfcP = 2, and the 
nonlinear Schrodinger equation follows from the many- 
particle dynamics as 



( 30 ) 
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FIG. 4. (Color online) Many-particle (blue [dark gray]) and 
mean-field (magenta [light gray]) energies in dependence on e 
for V = 1 and N = 30 particles. 


In this version the probability to find the system in the 
atomic mode is given by [xa|/ 2 , while the probability 
to find the system in the molecular mode is |xb|^. The 
variables sj are then defined as 


^ {X*aXb + X*bXa) , 

Sy = iXaXb - xlxa) , (31) 

= \ (iXal -2[X&|^) . 


Their dynamics are again given by equation (14). 


IV. MEAN-FIELD AND MANY-PARTICLE 
CORRESPONDENCE 

Let us first compare the spectral features of the mean- 
field and many-particle descriptions. For this purpose 
we define the mean-field energies as the values of the 
Hamiltonian function in the mean-field fixed points, 
that correspond to stationary solutions of the nonlinear 
Schrodinger equation. The resulting mean-field energies 
are plotted as a function of e for u = 1 in comparison 
with the many-particle energies for N = 30 particles in 
figure It can be clearly seen how the pattern of nar¬ 
row avoided crossings in the many-particle spectrum is 
closely following one of the mean-field energies. This en¬ 
ergy corresponds to the unstable fixed point at the tip of 
the teardrop that is associated to the all-molecular mode. 
At the critical values of e = ±-\/2 this becomes a stable 
elliptic fixed point, that is associated to the minimum 
and maximum eigenvalues of the many-particle system, 
respectively. The maximum and minimum mean-field en¬ 
ergies bound the many-particle spectrum, if the latter 
is renormalised by the semiclassical parameter 77 , as is 




FIG. 5. (Golor online) Mean-field (black) and many-particle 
dynamics for v = 1 and e = 1 (top) and e = 0 (bottom), and 
different particle numbers. The blue (dark gray), green (gray), 
and yellow (light gray) curves correspond to A = 20 , 100 , and 
500, respectively. The initial states are chosen as the ground 
states of TAa; in the top figure and —K^ in the bottom figure, 
and the initial conditions are marked by a star. The right 
figure on the top shows the dynamics of the ^-component for 
e = 1 , and the figure on the bottom right shows the dynamics 
of the z-component for e = 0 . 


expected in a typical quantum-classical correspondence. 
This is very promising for a semiclassical quantisation 
that we shall attempt in the next section. 

Let us now briefly turn to the correspondence between 
mean-field and many-particle dynamics. In figure we 
show several examples of mean-field trajectories (black 
lines) for different parameter values and initial condi¬ 
tions, in comparison to the corresponding many-particle 
dynamics. For comparison the initial many-particle state 
has been chosen as the ground state of a Hamiltonian of 
the type 

k = ak^ + bK, +cky, (32) 

whose expectation values of lie as close as possi¬ 

ble to the teardrop surface of the mean-field system. In 
figure 1 ^ we show the expectation values of k^ and k^, 
for this family of states for c = 0 and varying values of 
a and b, for different particle numbers. It can be seen 
how the mean-field teardrop is approached with increas¬ 
ing particle number. See m for a comparison between 
these states and a new type of coherent states proposed 
there. 

In figure it can be seen how the many-particle dy¬ 
namics closely follows the mean-field dynamics even for 
relatively small particle numbers in the vicinity of ellip¬ 
tic fixed points. Where the influence of the hyperbolic 
fixed point at the bottom of the teardrop is stronger, 
the deviations become larger and we observe the typical 
breakdown phenomenon of the many-particle dynamics. 
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FIG. 6. (Color online) Expectation values of and for 
the ground states of the Hamiltonian (321 where c = 0, b 


varies from —0.5 to 0.5 and a = ±|\/(1 — 26) (1 + 25)^, for 
different particle numbers N = 2 (blue dashed dotted line), 
A'’ = 4 (red dashed line), = 10 (magenta dotted line) and 
N = 100 (black solid line). 


FIG. 7. (Color online) Mean-field phase space portraits (top) 
and potential curves (bottom) for v = 1 and e = 0 (left) and 
e = 2 (right). The false colours in the phase space plots rep¬ 
resent the energy values. An orbit belonging to one selected 
energy each (marked by a black line in the bottom pictures) 
is highlighted in each of the phase space plots by a thick black 
line. 


For longer times (not depicted), revival phenomena occur 
where the revival time rapidly increases with the parti¬ 
cle number as expected m- The slow convergence of the 
many-particle dynamics towards the mean-field dynamics 
is particularly pronounced for those initial states whose 
mean-field dynamics asymptotically approaches the sad¬ 
dle point related to the all-molecular mode, as seen at the 
bottom in figure The breakdown happens very rapidly 
here and even for longer times, no revivals are observed. 
This behaviour has been investigated in some detail in 
[7], and is typical for quantum-classical correspondence 
in the neighbourhood of hyperbolic fixed points. 


V. SEMICLASSICAL QUANTISATION 

In the following we shall investigate how semiclassical 
techniques might be modified to obtain an approxima¬ 
tion for the many-particle eigenvalues from the mean- 
field dynamics. For this purpose it is convenient to use 
the approach by Braun [55] that has been developed for 
general three-term recurrence relations, since the matrix 
representation of our many-particle Hamiltonian ([^ in 
the eigenbasis of Kz is tri-diagonal, and thus the eigen¬ 
value equation defines a three-term recurrence relation. 
The same procedure has been successfully applied to 
the spectrum of the two-mode Bose-Hubbard model in 
[IIEIIS]. Since the spectrum of Kz is related to that 
of two coupled harmonic oscillators, we expect that the 
semiclassical quantisation should yield exact results for 
the spectrum of Kz for arbitrary particle numbers. For 
this purpose the quantisation procedure of Braun has to 
be slightly modified. We shall give a brief overview over 
this modified version in what follows. 


The main idea is to derive a WKB type approximation 
for the wave function in dependence on p, the match¬ 
ing conditions between the allowed and forbidden regions 
then yield a Bohr-Sommerfeld type quantisation condi¬ 
tion for the eigenvalues. The latter quantisation con¬ 
dition can in fact be motivated heuristically, by argu¬ 
ing that every many-particle state takes up on average 
a phase space volume of h, where the ground state and 
the highest exited state only take up the minimum uncer¬ 
tainty volume of |. However, the role of the semiclassical 
parameter h in the present context is taken over by the in¬ 
verse matrix size of the problem, that is, 77 = ("y -I- l) 
The quantisation condition then reads 

SipEn) = 2TTr]{n +^), (33) 

where S{E) is the mean-field phase-space area enclosed 
by the orbit of energy E. To calculate these phase space 
areas it is convenient to introduce potential curves for 
p. These potential curves U~^ and U~ are defined as the 
maximum and minimum functions of the Hamiltonian 
with respect to the angle dimension [55]. That is, we 
have 

C/± (p) = ep± ^ V(l-2p)(l + 2p)2, (34) 

which join at p = ±^, at the values ±|. Two examples 
of the potential curves together with the corresponding 
phase-space portraits are shown in figure The poten¬ 
tial curves bound the possible values of the energy of 
the mean-field system for a given value of p. On the 
other hand, for a given value of the energy, the motion 
is restricted to values of p where U~ (p) < E < C7+ (p), 
i.e. the classically allowed region. This region is bound 
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by the classical turning points p±, where E = U^{p±), 
which are found as the two roots of the polynomial 

2uV + “ (Y ^ “ T 

that fall into the interval [—|]. For any allowed value 
of the energies the polynomial has three real roots, two 
of which are located in the interval ^]. The third 
root, po, lies outside the physically relevant interval, to 
the left, that is, po ^ ^ P-- The fixed points of 

the dynamics correspond to the extrema of the potential 
curves and the tip at which both potential curves meet 
at p = — I. For the extremal values of the energy the two 
turning points coalesce. 

The phase-space area enclosed by an orbit with energy 
E can be written as the integral 

S{E) = j) q{p,E)dp, (36) 


where q{p) follows from the conservation of the energy as 


q{p) = arccos 


2 (g - ep) 
Y(l- 2 p)(l + 2 p )2 


(37) 


and care has to be taken that the area enclosed by the 
curve is calculated, rather than the area outside the 
curve. Depending on which of the potential curves which 
turning point lies on, the enclosed area in equation (36) 
can be evaluated as 


27r(p+-p_) - 2,§(F;), p±onC/_, 


S{E) = 


with 


\2n{\-pJ)-2~S{E), 

\2tt{\+p+)-2S{E) 

[ -2tt + 2S{E) 


P- on [/_, p+ on C7_|_, 
p_ on [/+, p+ on [/_, 
p± on U+, 

(38) 


rP+ 


SiE)= / q{p)dp. 


(39) 


In the case v = 0 this can be evaluated analytically, 
and we obtain the exact eigenvalues for arbitrary par¬ 
ticle numbers, as expected. For non-zero values of v 


the quantisation condition (33) can be straight forwardly 


solved numerically to obtain the semiclassical many- 
particle spectrum from the mean-field dynamics. In fig¬ 
ure we show the resulting semi-classical eigenvalues in 
dependence on e for v = 1 and two different particle 
numbers, in comparison to the numerically exact many- 
particle energies. We observe that even for iV = 4 par¬ 
ticles, corresponding to a relatively large value of 77 = |, 
the agreement is very good. As expected the quality of 
the semiclassical quantisation gets better with larger par¬ 
ticle number and with larger e, as in the limit v/e —>■ 0 
the quantisations yields exact results. 

The quantisation condition can be used to derive an 
analytic expression for the many-particle density of states 




FIG. 8. (Color online) Many-particle (solid blue lines) and 
semiclassical (dashed magenta lines) eigenvalues in depen¬ 
dence on e for v = 1 and two different particle numbers. The 
upper panel corresponds to A = 4 particles, the lower to 
N = 20. 


in the semiclassical limit of large particle numbers m 
[55] . Differentiating the quantisation condition (33) with 
respect to the energy yields 


dS'(A„) dn 


(40) 


That is, the many-particle density of states is given by 
the derivative of the mean-field phase space area with 
respect to the energy. This, on the other hand is given 
by the period of the mean-field orbit of the given energy 
[29) . That is, we have for the many-particle density of 
states 


dn _ T(E) 
dE ~ 2tt ’ 


(41) 


where T(E) denotes the period of the orbit with energy 


















E. This can be directly calculated from 


. , r+ dp fP+ fdH\ 

= (w) 

P+ dp 


= 2 


-1 


dp 


'p 


^{U+ -E){E-U-y 


(42) 


From the explicit form of the potential curves we hnd the 
period in terms of an elliptic integral as 


T{E) = ^ 


2 V 2 


de 


Jo 




P+-PO 


(43) 


P+ -P- 


’Wp+ - Po \P+-Po 


where po is the third root of the polynomial (35), and 


K denotes the complete elliptic integral of the first kind. 
The period diverges at the classical turning point, which 
can be seen as follows. At the value A = — |, correspond¬ 
ing to the orbit passing through the tip of the teardrop, 
the polynomial (35) can be explicitly factorised as 


PiE = -|) = (p + ^yi2v‘^p - 


+ a- 


(44) 


That is, in the subcritical case |e| < V^lvl, where the 
fixed point at the tip corresponds to a saddle point, we 
have Po = P- = —y and p+ = ^ — that is, we have 
T{E) oc K{1), which diverges. For supercritical values, 
on the other hand, where the fixed point corresponds to 
the minimum of the energy, we have p- = p+ = ~ 5 > and 
1 2 1 

« ^ 1 period has the finite value 


Po= 2 
of^ 




(p < - 3 . Thus, 
in this case. 


Figure depicts the normalised histogram of many- 
particle eigenvalues for a large particle number of = 
10000 in comparison to the mean-field periods (divided 
by 27r), which are given by the analytical expression 
(43). Note that the mean-field energies are rescaled with 
respect to the many-particle energies by rj. An excel¬ 
lent agreement between the analytical expression and the 
many-particle density of states is observed. In particu¬ 
lar the accumulation of states around the classical saddle 
point for values of e below the critical value is nicely re¬ 
covered in the many-particle histogram. Since the period 
of the orbit through the saddle point is infinite, the accu¬ 
mulation of many-particle eigenstates at this point leads 
to an actual divergence in the limit A^ —>■ 00 , which can 
be connected to a quantum phase transition [5]. 

We can also obtain an approximation for the compo¬ 
nents of the eigenvectors in the standard basis of K^, via 
a simple WKB ansatz. In the classically allowed region 
(between the two turning points) we make the ansatz 
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FIG. 9. (Color online) Density of states ^ for the many- 
particle system in comparison with the mean-held periods for 
V — 1 and different values of e. The many-particle density 
of states is approximated by the normalised histogram of the 
energies, for N = 10000 particles. The mean-held result is 
depicted by the solid magenta line. The values of e are e = 
0, 1, 2, 5 from top left to bottom right. 


where Wd denotes the classical probability distribution 
-1 


\Wcl {P)\ = 


1 


dH 

dq 


1 


2T^f{l-2p)il + 2pf-{E-epy 
and the action S{p) is given by 


.,(46) 


S{p) = 


with 


'■k{p-P_) - S{p), 

'^ih+P) ~ 

—TT - 1 - S{p) 


p± on U-, 

P- on U-, p+ on U+, 
P- on [/+, p_|_ on [/_, 
p± on C/+, 

(47) 


S{p)= / q{p)dp. 


(48) 


In the forbidden region the WKB ansatz reduces to the 
single exponential decaying solution 


IV'ci {p)\^ = ^ 


2i 

ujci {p) exp (- S{p) 

V 


(49) 


with 


T /p q (,p) dp, p<p-,p- onU±, 
T q (P) dp, p> p+,p+ onU±. 


S{p) = 

The matching conditions at the boundary then impose 
the quantisation condition (33) and the absolute value of 
the WKB wave function in the classically allowed region 
becomes 


i’ip) = VwdXp) -I- A_e , (45) 


\ipci {p)\^ = ‘2\wciip)\cos^ (jl ^'S'(p)-0- (51) 
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FIG. 10. (Color online) Exact n-th many-particle eigenvec¬ 
tors (green circles) with the WKB wave function (red) for 
N = 40, e = 0.5 and v = 1 for n = 1 (top), n = 3 (middle) 
and n = 10 (bottom). The red crosses indicate the semiclas- 
sical wave function ( |52[ ). The solid red line, indicating the 
continuous semiclassical wave function according to ( |51[ ), is 
added to guide the eye. 


While this expression appears to depend on the contin¬ 
uous variable p, for finite values oi t] p can only take 
on discrete values, due to the periodicity of its conjugate 
variable q, to which it can be related via a discrete Fourier 
transform. We thus have the semiclassical approximation 

l^'n) = i’ci{v'm)\m), (52) 

m=-f :l:f 

where \m) denotes the eigenvectors of Kz, that is, the 
states with 2m + ^ atoms and ^ — m molecules. Fig¬ 
ure [^depicts examples of the exact many-particle eigen¬ 
vectors in comparison to the semiclassical approximation 
(normalised to fit the central maximum). As expected, 
the semiclassical approximation breaks down in the vicin¬ 
ity of the turning points, but approximates the many- 
particle wave functions well for other values. 


VI. SUMMARY AND OUTLOOK 

While the role of many-particle effects in cold atom 
systems is crucial, in large realistic systems the mean- 
field approximation is often all that is accessible. Thus, 
the possibility to recover many-particle features from the 
mean-field description is an important addition to the 
theoretical toolbox for cold atoms. Here we have demon¬ 
strated that semiclassical methods can be modified to 
deduce many-particle properties for atom-molecule con¬ 
version systems from the mean-field dynamics alone. We 
have considered the eigenvalues and eigenvectors here; an 
extension to dynamical properties via semiclassical prop¬ 
agators is an important topic for future investigations. 
A non-trivial issue is the generalisation to more realistic 
models with many modes for both atoms and molecules. 
Progress has been made in this direction for cold atomic 
systems without atom-molecule conversion in [HH]. The 
combination with the results obtained here suggests that 
this goal is not out of reach. An obstacle for the ap¬ 
plication of some semiclassical techniques is the absence 
of a well-defined set of condensed states that coincide 
with coherent states of the deformed SU{M) algebra for 
atom-molecule conversion systems. There have been pro¬ 
posals for coherent states in misni, however they do not 
fulfil all the properties that one relies on in the case of 
cold atoms in M mode systems when employing SU{M) 
coherent states m- The connection between coherent 
states and projective manifolds [32] for SU{M) systems 
might lead a way forward. 
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